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We calculate the pion's valence-quark momentum-fraction probability distribution using a Dyson- 
Schwinger equation model. Valence-quarks with an active mass of 0.30 GeV carry 71% of the 
pion's momentum at a resolving scale go ~ 0.54 GeV — 1/(0.37 fm). The shape of the calculated 
distribution is characteristic of a strongly bound system and, evolved from go to g = 2 GeV, it 
yields first, second and third moments in agreement with lattice and phenomenological estimates, 
and valence-quarks carrying 49% of the pion's momentum. However, pointwise there is a discrepancy 
between our calculated distribution and that hitherto inferred from parametrisations of extant pion- 
nucleon Drell-Yan data. 

Pacs Numbers: 13.60. Hb, 14.40.Aq, 12.38.Lg, 24.85. -fp 



I. INTRODUCTION 

The cross section for deep inelastic lepton-hadron scat- 
tering can be interpreted in terms of the momentum- 
fraction probability distributions of quarks and giuons 
(partons) in the hadronic target, and since the pion is a 
two-body bound state with only u- and c?-valence-quarks 
it is, in some respects, the simplest hadron and therefore 
represents the least complicated, nontrivial system for 
which these distribution functions can be calculated. In 
a theorist's ideal world, stable pion targets would then fa- 
cilitate a comparison between these calculations and the 
parton distribution functions measured in deep inelastic 
scattering. However, pion targets are not abundant, and 
therefore these functions have primarily been inferred 
from Drell-Yan and direct photon production ||] 

in pion-nucleon and pion-nucleus collisions: an approach 
that can only be successful if the nucleon's parton distri- 
butions are well known |4| . More recently, semi- inclusive 
ep eNX reactions, N = p,n, have also been em- 
ployed , a method which assumes that small- virtuality 
pions dominate leading nucleon production. These ex- 
periments are challenging but high-statistics data exist. 

The theoretical situation is arguably poorer: the sys- 
tematic tool provided by perturbative QCD does not ad- 
mit the calculation of the distribution functions but only 
of their evolution from one large spacelike-g^ to another; 
and numerical simulations of lattice-QCD are currently 
restricted to the quenched approximation and only yield 
moments of the distributions, not the distributions them- 
selves [||. Furthermore, the fact that the pion is both 
a bound state and the Goldstone mode associated with 
dynamical chiral symmetry breaking complicates the cal- 
culation of the pion's distribution functions and places 
additional constraints on any framework applied to the 
task. Considerations of chiral symmetry have led some 
[0-^ to adopt the Nambu-Jona-Lasinio model as a ba- 
sis for their calculations but a number of the difficulties 
with this approach, among them a marked sensitivity to 




the regularisation procedure in this non-renormalisable 
model, are emphasised and discussed in Refs. 
stituent quark models have also been employed | IC 
difficulties encountered in such studies considered in Ref. 
pd| ; and so has an instanton-liquid model p ^ . 

The Dyson-Schwinger equations (DSEs) |p,3| provide 
an approach well-suited to the calculation of pion ob- 
servables. Since a chiral symmetry preserving trunca- 
tion scheme exists they easily capture the dichoto- 
mous bound-state/Goldstone-mode character of the pion 
[|l5| , ^ . Furthermore, because perturbation theory is re- 
covered in the weak coupling limit, they combine; e.g., a 
description of low-energy tt-tt scattering |17| with a cal- 
culation of the electromagnetic pion form factor, F.^{q^), 
that yields (l8|: the 1/g^-behaviour expected from per- 
turbative analyses at large spacelike-g^ and a calculated 
evolution to the p-meson pole in the timelike region [ p^ . 
These and other features of contemporary applications 
are described in Refs. 

Herein we employ a phenomenological DSE model, 
used previously |22p3] in a description of a wide range 
of electron- nucleon and meson- nucleon form factors, in 
a calculation of the pion's valence-quark distribution. 
The framework is Poincare covariant and exhibits, in 
a simplified form, the features of renormalisability and 
asymptotic freedom characteristic of QCD. Another im- 
portant aspect of our study is a description of the pion 
as a bound state of a dressed-quark and -antiquark 
through its momentum- rfependent Bethe-Salpeter am- 
plitude, with the dressed-quark propagator exhibiting 
that momentum-dependence characteristic of DSE stud- 
ies and recently confirmed in simulations [ p4[ of lattice- 
QCD. 

Our article is organised thus: Sec. H recapitulates the 
calculation of the hadronic tensor for lepton-pion scat- 
tering; Sec. HI describes the DSE elements necessary in 
the calculation of the distribution functions - dressed- 
propagators, Bethe-Salpeter amplitudes, etc.; Sec. IV re- 
ports and discusses the results; and Sec. V is an epilogue. 
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FIG. 1. "Handbag" contributions to the virtual pho- 
ton-pion forward Compton scattering amplitude, which are 
the only impulse approximation diagrams that survive in 
the deep inelastic Bjorken limit, Eq. (^). tt, dashed-line; 7, 
wavy-line; S, internal solid-line, dressed-quark propagator, 
Eq. (|o|). The filled circles represent the pion's Bethe-Salpeter 
amplitude, F-^ in Eq. (jssj), and the dressed-quark-photon ver- 
tex, F^ in Eq. (p7|), depending on which external line they 
begin/end. 



II. LEPTON-PION SCATTERING 

Deep inelastic scattering from a pion target can be 
studied via the diagrams in Fig. ^ which describe virtual 
photon-pion forward Compton scattering and are p5[ | the 
only impulse approximation contributions to survive in 
the Bjorken limit:^ 



CX3 , P ■ q 



-00 but X 



2P-q 



fixed. (1) 



The upper diagram represents the renormalised matrix 
element 



f 



(2) 



■ T_r7r(fcr; ~P) S{kt) ieQTv{kt, k) 

(27r)* 

X S{k) ieQr^{k, kt) S{kt) r+r^(fcr; P) S{k,) , 



'^We use a Euclidean metric convention, in which a ■ b :— 
ifj. bv Sf^i, :— 'Yli=i '^i ' that a spacelike vector, Q^, has 



Q > and an on-shell pion is described by P 



Our 



Dirac matrices are Hermitian and are defined by the algebra 



where: 
and 



Ttt (£; P) is the pion's Bethe-Salpeter amplitude 



f^{t-P) = C^T^{-t-P)'^C , 



(3) 



with T± — i(Ti ± iT2), C — 7274, the charge con- 
jugation matrix, and denoting matrix transpose; 
S{t) — diag[S'„(.^), S'£i(^)] is the dressed-quark propaga- 
tor and we assume Su = Sd = S throughout (isospin 
symmetry of the strong interaction); r^(£i,^2) is the 
dressed-quark-photon vertex, with Q = diag(2/3, — 1/3) 
the quark-charge matrix; kr = k — q — P/2, kt = k — q, 
ks = k — q — P; and the trace is over colour, flavour and 
Dirac indices. (In Eq. H, since the renormalised matrix 
element is finite in our DSE model, as in QCD, we have 
not made explicit the usual translationally invariant reg- 
ularisation scheme [^ .) The matrix element represented 
by the lower diagram is the crossing partner of Eq. (j^) 
and is obvious by analogy. 

The hadronic tensor relevant to inclusive deep inelastic 
lepton-pion scattering can be obtained from the forward 
Compton process via the optical theorem: 

^) = ^I-^ P) + T-^T, P)] , (4) 

and because of current conservation it may be expressed 
in terms of only two invariant structure functions: 

W,,{q- P) ^W,{q\P- q) t,, - ^^^^'"^ ' , (5) 

with tf,i, = - q^,q^/q^ and Pf, ^ Pf, - q^, P ■ q/ q^ ■ 

To proceed it is useful to express the dressed-quark 
propagator as 



s{k)^ ^{e)^^{k), 

A(fc2) ^ 1/[A:2 + M2(fc2)], 
N{k) ^ Z{k^) [-^7 • k + M{k^)\ , 



(6) 
(7) 
(8) 



where Z{k'^) is the dressed-quark wave-function renor- 
malisation and M{k'^) is the dressed-quark mass func- 
tion. Using Eq. (^) and evaluating the colour and flavour 
traces, Eq. (||) assumes the form 



r+(9,P) 
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J^Aie)A{k',)T+{k,q,P), 



(9) 



where 



T+ =. tTD [f ,(fcr; -P) Sih) iT,ih,k) 
X Mik) iT^ik, kt) Sih) r,(fcr; P)M{ks)] 



(10) 



and the trace is now, of course, only over Dirac indices. In 
this form it is easy to adapt the procedure of Ref. and 
isolate the pinch singularities that yield the imaginary 
part in Eq. (^). 

Introducing the integration variable transformations: 
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k = K + yq + zP, K ■ q = = K ■ P , 



(11) 



which has Jacobian J = P ■ q in the Bjorken hmit, Eq. 
(0), and subsequently 



J/=l 



:P2 



2P-q 2P-q' 



one finds 



k"^ = 2P-q{z-x) + 0{{P-q)°) 



(12) 



(13) 



so that on the domain relevant to inclusive deep inelastic 
scattering 



1 



1 



2P-qz~x 2P-q 



S{z — x) 



(14) 



and the contribution of r+ to Wp,y can be written 



1 1 



2TTAP-q 9 



(2 



^ 1 dwd?K lm[A(fc2)T+(fc,g,P) 



(15) 



This shows that in the Bjorken limit the struck parton 
carries a fraction x of the pion's momentum. 

The integrand can now be simplified further by using 
Eqs. (pi]) and (|l^) to express 



(16) 



s:^ki^{x- 1) [w - + K^/ix - 1)] 
and hence 



A(.) 
so that 



1 



1 



\-x w-P-^-[K-^ + M-^{s)]l(l - x) 



(17) 



1 



1 



P -q l-x 



1 



<fK -[T+ik,q,P) 



w—w{k) 



(18) 
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with w{k) — P^ + (k^ + M^)/(l — x), where the valence- 
quark mass is obtained as the solution of M = M{—M'^). 

The momentum of the remaining dressed-quark line 
also simplifies: 



ji := kf = 



(19) 



which illustrates that the integrand depends only on 
as an integration variable. Another shift of variables is 
now useful: J (Pk = tt J diP , 



= (1 - x) (^ - ^min) , 

Mmin = a;[P'+MV(l-a;)], 



(20) 
(21) 



and hence 



1 



(47r)2 



1 



4P-g 



[T+{k,q.P)\ 



(22) 



with w{^) = M + Aimin(l — x)/x. 

It is apparent in Eq. (^2|) that the only remaining func- 
tional dependence that is not obviously expressed solely 
in terms of x is that arising from the nontrivial aspects 
of the Dirac trace in Eq. (M). Nevertheless, direct cal- 
culation; e.g., shows that in the Bjorken limit 

[T^uik, q, P)] ir:;^"^ ^P-q {a{x) t,^ + . . .) (23) 
so that one has 



W+{q;P)^F+{x)t^.- 



and the Callan-Gross relation: 



P-q 



F+{x) = 2xF+{x) . 
At this point it is obvious that 

Pi 2 (a;) ^ as a; ^ 1 



(24) 



(25) 



(26) 



because of the contraction of the integration domain: 

/^iniii ^ as X 1. The analysis of W^^, follows a 

similar pattern and yields the same results, and Fi 2 — 

F+ 
1,2 



1,2- 



Equations (^), (10), ( p^ and ( p5| ) provide a model- 
independent starting point for calculating the valence- 
quark distribution functions using model input for the 
internal elements represented in Fig. |l]; i.e., the dressed- 
quark propagator, pion Bethe-Salpeter amplitude and 
dressed-quark-photon vertex. These are valence-quark 
distributions because, although sea-quarks are implicitly 
contained in the dressing of the propagators and calcula- 
tion of the dressed vertices, the "handbag" impulse ap- 
proximation diagrams in Fig. |^ only admit a coupling of 
the photon to the propagator of the dressed-quark con- 
stituent. The internal structure of the dressed-quark is 
not resolved and therefore the calculation yields the dis- 
tribution at a scale qo characteristic of the resolution: 
qo is an a priori undetermined parameter in calculations 
such as ours, although we anticipate 0.3 < qo ^ 1-0 GeV, 
with the lower bound set by the Euclidean constituent- 
quark mass and the upper by the onset of the pertur- 
bative domain. A sea-quark distribution is generated 
via the renormalisation group (evolution) equations when 
the valence distribution is evolved to that g^-scale appro- 
priate to a given experiment. To generate explicit sea- 
quark contributions at the scale requires going beyond 
the impulse (handbag) approximation; e.g., incorporat- 
ing photon couplings to the intermediate-state quark- 
meson-loops that can appear as a dressing of the quark 
propagator: 7r+ = u d — + (uss) d = (K^s) d u d = tt^ , 
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with jK'^s K'^s'j, etc. Such intermediate states arise 
as vertex corrections in the quark-DSE. 
With these observations in mind then 

P2{^:Qo) = ^xu^(x;qo) = |x[M(x;qo) - u{x;qo)] , 
F2{x;qo) = ^xdy{x;qo) = lx[d{x;qo) - d{x;qo)] , 

(27) 

and it is straightforward to demonstrate algebraicaUy 
that 



60 bi 



The calculations are required to yield 



1 /•! _ 

dxUv{x;qo) = l= / dxdv{x;qo): 
^0 



(28) 



(29) 



i.e., to ensure that the 7r+ contains one, and only one, 
M-valence-quark and one d-valence-quark. 

We remark that (x) — d'^ [x) in the fj-parity sym- 
metric limit. However, in this limit lo — s- tttt is forbidden 
and hence the scale of tj-parity symmetry violation in na- 
ture is characterised by the ratio ^ uj^tttt p^ir-r; — 
0.1%. This bound on any difference between the pion's 
quark distribution functions is consistent with the model 
estimate in Ref. E7||. 



III. MODEL ELEMENTS 

To complete our calculation the internal lines and irre- 
ducible vertices appearing in Fig. ^ must be specified. In 
general these quantities can be obtained by solving the 
quark DSE and the appropriate inhomogeneous Bethe- 
Salpeter equations | |20[ |. However, the study of an exten- 
sive range of low- and high-energy light- and heavy-quark 
phenomena has yielded efficacious algebraic parametrisa- 
tions p^ , p3| , p8| -^ , and we employ those herein. 

The dressed-quark propagator is 



= [z7-pA(p2) + S(p2)]-i ^ 

as{x) = 2mT[2{x + m^)) 

^'T(hxx) T[h-ix) [60 + h2T[ex)\ 
1 



ay(x) = 



[1- J^(2(a; + m2))] 



(30) 

(31) 
(32) 

(33) 



with T{y) = (1 — c^y)/y, x — p^/A^, rfi — mlX, as{x) — 
^<^s{p^) and <fv{x) — A^cry(p^). The mass-scale, A = 
0.566 GeV, and dimensionless parameter values:n 



^ e — 10^* in Eq. (^) acts only to decouple the large- 
and intermediate-p^ domains. The study used Landau gauge 
because it is a fixed point of the QCD renormalisation group 
and Z2 ~ 1, even nonperturbatively [Q. 



0.00897 0.131 2.90 0.603 0.185 



(34) 



were fixed in a least-squares fit to light-meson observ- 
ables and the dimensionless u-current-quark mass 
corresponds to 



(35) 



This algebraic parametrisation combines the effects of 
confinement and dynamical chiral symmetry breaking 
with free-particle behaviour at large spacelike |2C| ] 
and, as illustrated explicitly in Ref. |^2|, its qualitative 
features have recently been confirmed in simulations of 
lattice-QCD ||. 

Since the quark described by Eqs. (|30|)-(|3^) is con- 
fined there is no solution of the mass-shell equation: 
,s + Af^(s) = 0. However, in calculations of observables, 
the Euclidean constituent-quark mass, M^, defined as 
the solution of s = M^(s), provides a realistic estimate 
of the quark's active quasi-particle mass. In the present 
example. 



0.33 GeV. 



(36) 



The large value of the ratio: /m^'^'^^ fa 70 is one 
signal of the effect that the dynamical chiral symme- 
try breaking mechanism has on light-quark propagation 
characteristics. Another is the value of the chiral-limit 
vacuum quark condensate and one advantage of using the 
algebraic parametrisation is that it yields the following 
simple expression p8|: 



47r2 bi &3 ' 



(37) 



which makes plain the relation between the condensate 
and the chiral-limit coefficient of the l/p'*-term in a sip) 
@- 

The general form of the 7r-meson Bethe-Salpeter am- 
plitude is 

r.(fc;Q)-75 iE^{k;Q)+j-QF^{k;Q) 



J ■ kk ■ Q G-^{k;Q) + a^u k^Qy H^^ik; Q) 



(38) 

where the behaviour of the invariant functions is con- 
strained to a large extent by the axial-vector Ward- 
Takahashi identity jl^Jl^. That identity and numerical 
studies m& support the parametrisation 



EAk;Q) = ^B4k^) 



(39) 



where B^^ is obtained from Eqs. (|30[)-(p3|), evaluated us- 
ing rh — and 



60 ^ 6n = 0.19 



(40) 
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with the other parameters unchanged, and: 



F^ik;Q)^E^ik;Q)/{nOU) 
G^{k-Q)^2F^{k-Q)/[k^ 



(41) 
(42) 



M\]Y = IOAqcd; and HTr{k;Q) = 0. The amphtude is 
normahsed canonicaUy and consistent with the impulse 
approximation {k± — k ± Q/2): 



d^k 

(2^ 
+tr|?' 



tr^ 



f-(fc;-Q)S^r,(fc;Q)^(fc-) 



dQn 



T^{k;^Q)S{k+)r^ik;Q) 



dS{k^ 



(43) 



with the trace over flavour indices omitted, and this fixes 
Nt^. The often- neglected pseudovector elements: i^^, G,r, 
are crucial to recovering the 1/q^ behaviour of the 
electromagnetic pion form factor at large spacelike-g^ ex- 
hibited in perturbative QCD analyses. However, models 
efficacious for low pion-energy can be constructed with- 
out them |l||l,|0|j3l]j3|]. 

In Ref. ||15|| a pseudoscalar meson mass formula 
was derived that, in the limit of small current-quark 
masses, reproduces what is commonly known as the Gell- 
Mann-Oakes-Renner relation, and also has an important 
corollary applicable to mesons containing heavy-quarks 
1^,^. In the case of the pion it gives 



r2 2 



(44) 



with the parametrisations yielding an algebraic expres- 
sion for the "in-pion" condensate: 



(g-<z)^ = ln(CVA|cD) 



(45) 
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and where the leptonic decay constant is obtained from 



fyrQu NctTD 



{2n)' 



■-f5-f^S{k+)TAk;Q)S{k^). (46) 



It remains to specify the dressed-quark-photon vertex. 
The manner in which an Abelian gauge boson couples to 
a dresscd-fcrmion has been much studied and a range of 
qualitative constraints have been elucidated This 
research supports an Ansatz [p6[ 



j:F{el,el) = i[F{el) + F{il)]. 



F{ll)-F{il) 



(47) 
(48) 

(49) 



where F = A,B are the scalar functions in Eq. (30), 
which preserves many of the constraints, important 
among them the vector Ward-Takahashi identity, and 



is expressed solely in terms of the dressed-quark prop- 
agator. In concert with the algebraic parametrisa- 
tions described above, it has been widely used in the 
study of electromagnetic processes and is phenomeno- 
logically efficacious; e.g., providing for the parameter- 
independent realisation of "anomalous" photon-hadron 
couphngs mm^]. For these reasons, we employ Eq. ( |47| ) 
herein. Nevertheless, significant progress has recently 
been made with the direct calculation of 

in DSE models of QCD and those studies will provide 
the basis for improved Ansdtze in the future. 

The elements described herein yield the following cal- 
culated values for a selected, illustrative range of light- 
hadron observables: 





Calc. 


Obs. 




/.(GeV) 


0.090 


0.092 


26 




0.139 


0.138 


26 


(-(g-g)0)i/3 


0.221 


0.236 


40 




0.250 


0.245 


16 




0.56 


0.663 


41 




0.78 


0.87 




^r2(fn?) H 


(0.33)2 (0.34)2 



(50) 

with r^r calculated in impulse approximation (see, e.g., 
Refs. jl^,^) and the condensates evaluated at C = 
1 GeV using a 3-flavour value of Aqcd — 0.242 GeV. 



IV. CALCULATED DISTRIBUTION FUNCTIONS 

We can now proceed with the evaluation of the one- 
dimensional integral that yields Uy{x) via Eqs. ( p2[ ) and 
(p7|). That calculation requires a determination of the 
valence-quark mass, M. The parametrisation of the 
dressed-quark propagator is confining and does not admit 
a solution of M = M{—M'^). However, in inclusive deep 
inelastic scattering, confinement is recovered through in- 
coherent hadronisation after the dissolution of the bound 
state and we introduce that aspect herein by adopting a 
quasiparticle representation: A(fc2) — l/[k'^ + M^] in 
Eq. (^, with M determined by requiring that Eq. ( p9[ ) is 
satisfied. It is valence-quarks of this mass that populate 
the pion. This is an internally consistent prescription if 
M « M^. 

Our calculated form of Uy{x; go) is depicted in Fig. ||. It 
vanishes at a; = 1 , in accordance with the kinematic con- 
straint expressed in Eq. (|2^) , and corresponds to a finite 
value of Fi {x = 0), which is a signal of the absence of sea- 
quark contamination. Furthermore Uy{x = 0; go) 7^ Oi 
which is as it should be since there is no constraint that 
requires it to vanish at this point. Unsurprisingly for a 
light bound state of heavy constituents, the shape of the 
distribution is characteristic of a strongly bound system 
p3[ : cf. for a weakly bound system Uy{x) ~ 5{x — i). 
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X 

FIG. 2. Solid line: Uv{x;qo) calculated using the DSE 
elements and parameters described in Sec. The va- 

lence-quark mass is M = 0.30 GeV and the resolving scale 
qo = 0.54 GeV = 1/(0.37 fm) is fixed as described in con- 
nection with Eq. (p4). Dashed line: Uy{x;qo) calculated in 
the absence of the pseudovector components of the pion's 
Bethe-Salpeter amplitude; i.e., _F = = G in Eq. in- 
stead of Eqs. ( ^ ) and (^2|). Dot-dashed line: distribution 
calculated with — > 0.1 m^, M = 0.36 GeV. 

The area under each of the curves in Fig. ^ is one and 
requiring that for Uy(x;qQ) yields a calculated valence- 
quark mass of 

Af = 0.30 GeV, (51) 

which is within 10% of this model's Euclidean constitu- 
ent-quark mass, Eq. (|3^). The value of M affects the 
position of the peak in xuy{x): increasing M shifts the 
peak to lower x. 

The average momentum-fraction carried by the 
valence-quarks at this resolving scale is 




dxx[u.v{x; go) + d.v{x; qo)] = 0.71 , (52) 



with the remainder carried by the gluons that effect the 
binding of the pion bound state, which are invisible to the 
electromagnetic probe. The second and third moments 
of the distribution are 

=0.18, {x^)q„=0.10. (53) 

To determine the resolving scale, qo, we employ the 3- 
flavour, leading-order, nonsinglet renormalisation group 
(evolution) equations (see, e.g., Ref. |^^) to evolve the 
distribution in Fig. ^ up to q — 2 GeV, and require 
agreement between the first and second moments of our 
evolved distribution and those calculated from the phe- 
nomenological fits of Ref. jj] . With 

go = 0.54 GeV = 1/(0.37 fm) (54) 



X 

FIG. 3. Dashed line: xUv{x;qo); Short-dashed line: fit 
of Eqs. (H) and (|^; Solid line: the evolved distribution, 
xuv{x;q — 2 GeV); Dotted line: xuv{x;q — 4.05 GeV), 
evolved with a 4-flavour value of Aqcd ~ 0.204 GeV; and 
Dot-dashed line: the phenomenological fit of Ref. Q. The 
data are from Ref. obtained with an invariant /i^/i~-mass 
> 4.05 GeV and inferred from the differential pion-nucleon 
Drell-Yan cross section using simple distribution parametri- 
sations of the type in Eq. (BSl) . This data was part of the set 
employed in the fit of Ref. The fits of Eqs. (||) and (|^) 
are not depicted: they are almost indistinguishable from our 
calculated curves. 



we obtain at g = 2 GeV 











Gale. 


0.24 


0.098 


0.049 


Fit [4] 


0.24 ±0.01 


0.10 ±0.01 


0.058 ±0.004 


Latt. § 


0.27 ±0.01 


0.11 ±0.3 


0.048 ±0.020 



(55) 



at which point the valence-quarks carry a momentum- 
fraction of 0.49. At qo, as/(27r) = 0.13: this is the com- 
bination that appears in the evolution equation. (NB. At 
0.5 go, a,/(27r) = 0.68.) 

The original and evolved distributions are depicted 
in Fig. ^. The evident accentuation via evolution of 
the convex-up behaviour of the distribution near x ~ 1 
is characteristic of the renormalisation group equations, 
which populate the sea-quark distribution at small- a; at 
the expense of large- a; valence- quarks. 

A fit to xuy(x;q), adequate for the rapid estimation 
of moments, is provided by the simple data fitting form 
employed in Refs. H,^: 

xur-'ix; q) = (1 - x f , (56) 

r(a)r(i + /3) 

with in our case 

g (GeV) I 0.57 2.0 4.05 

a 1.34 0.92 0.84 (57) 
P 1.31 1.80 1.98. 
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The moments of u™°™-{x\ q) are given by 



V. EPILOGUE 



We emphasise, as is evident in Fig. ||, that Eq. (|5^) is 
not a good pointwise approximation to our calculated 
form of xUv{x\qQ) and, furthermore, Eq. ( p6[ ) divided- 
by X provides a very poor pointwise approximation to 
Uy(x] go) for x < 0.5. 

An alternative parametrisation has also been employed 
in fitting data: 

X = x^Ml - (1 ~euV^ + -fu x) , (59) 

with Au fixed by Eq. (p9[). This, with the parameter 
values: 



q (GeV) 


0.57 


2.0 


4.05 


Au 


11.24 


4.25 


3.56 


Vi 


1.43 


0.97 


0.89 


m 


1.90 


2.43 


2.61 




2.44 


1.82 


1.62 


In 


2.54 


2.46 


2.30 , 



(60) 



provides a pointwise accurate interpolation of the calcu- 
lated forms of X {x) depicted in Fig. ^: on the scale of 
this figure, the fit and our calculation are barely distin- 
guishable. Furthermore, for x > 0.2, M^"(x;qo) obtained 
from Eq. ( |5^ ) even provides a good pointwise approxima- 
tion to Uy{x; go) in Fig. ^ The difference between Eqs. 
( p6| ) and ( |59| ) is the improved capacity of the latter to ac- 
commodate convexity in the parton distributions, which 
is a signal feature of our calculation. 

Particular regularisations of the Nambu-Jona-Lasinio 
model yield a distribution with the functional form 
u^{x; Qq'^^) = 0{x) 9{l—x), which corresponds to valence- 
quarks carrying each and every fraction of the pion's mo- 
mentum with equal probability. That result is an arte- 
fact, arising from the representation of the pion bound 
state by a moTaentum- independent Bethe-Salpeter am- 
plitude; i.e., from representing the pion as a point- 
particle, which is a necessary consequence of the model's 
momentum-independent interaction. In this case, be- 
ginning at g^-'L = 0.35, for which a/(27r) = 0.31, 
X u^{x; Qq'^^) evolved to q = 2 GeV is pointwise very well 
described by Eq. (|^) with a = 0.67, (3 = 1.13, cf. the 
values inRef. g: a = 0.64 ± 0.03, /? = 1.15 ± 0.02. 

We infer from this result and the discussion above that 
the fitting form in Eq. ( |5^ ) is inadequate: it is contin- 
uously connected to a distribution that lacks dynamical 
content and is unable to represent that structure in the 
distributions which characterises dynamics. These ob- 
servations yield insight into the efficacy of the updated 
fitting form @ in Eq. @). 



We calculated the valence-quark distribution in the 
pion using a Dyson-Schwinger equation (DSE) model 
that provides a good description of a wide range of 
hadron observables: Figs. ^ and ^ summarise our new 
results. The DSE model describes valence-quarks with 
an active mass M — 0.30 GeV at a resolving scale 
qo = 0.54 GeV = 1/(0.37 fm), and with this value of go 
the evolution to values of the momentum scale relevant 
to contemporary experiments yields a distribution whose 
low moments agree with the values obtained in lattice 
simulations and from a phenomenological fit. 

There is a pointwise difference between our calculated 
distribution and the form used hitherto in parametris- 
ing pion data, as evident in Fig. ^. That discrepancy can 
also be observed in related covariant calculations ^,^,|l2| . 
From the information currently available, its origin does 
not appear to lie in model details, and we note that any 
assumed fitting form with little or no convexity in the 
vicinity of a; = 1, which is adapted to a body of data 
at a given scale, q^, will necessarily become convex- up 
under evolution to an higher scale, q^ > q^. Part or all 
of this discrepancy may therefore be attributable to the 
restricted function space used thus far in parametrising 
pion data. That possibility is supported by the capac- 
ity of an updated fitting form, used for the proton, to 
accommodate the structure evident in our calculations. 
However, we cannot rule out the possibility that the dis- 
crepancy may point to an as yet overlooked shortcoming 
in the application of models to the calculation of distri- 
bution functions. 

Our calculation is a prototype. It can be improved; 
e.g., by using direct numerical solutions of the quark- 
DSE and meson Bethe-Salpeter equations, as in recent 
calculations of the light-meson electromagnetic form fac- 
tors p^ , and/or incorporating sea-quark contributions 
at the soft-scale, q^. It can also be extended to more 
robust "targets," such as the nucleon, using models like 
those of Refs. p3|,i5|. 
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